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ABSTRACT 

Improved collisions strengths for the mid-infrared and optical transitions in Nev are 
presented. Breit-Pauli R-Matrix calculations for electron impact excitation are carried 
out with fully resolved near-threshold resonances at very low energies. In particular, 
the fine structure lines at 14 /im and 24 /ini due to transitions among the ground state 
levels ls'^2s'^2p^ ^^0,1,2, and the optical/near-UV lines at 2973, 3346 and 3426 A tran- 
sitions among the "^^0,1,2,^ D2,^ Sq levels are described. Maxwellian averaged collision 
strengths are tabulated for all forbidden transistion within the ground configuration. 
Significant differences are found in the low temperature range Te< 10000 K for both 
the FIR and the opitcal transitions compared to previous results. An analysis of the 
14/24 /im ratio in low-energy-density (LED) plasma conditions reveals considerable 
variation; the effective rate coefficient may be dominated by the very low-energy be- 
haviour rather than the maxwellian averaged collision strengths. Computed values 
suggest a possible solution to the anomalous mid-IR ratios found to be lower than 
theoretical limits observed from planetary nebulae and Seyfert galaxies. While such 
LED conditions may be present in infrared sources, they might be inconsistent with 
photoionization equilibrium models. 

Key words: Gaseous Nebulae - Optical Spectra: H 11 Regions - Line Ratios: Atomic 
Processes - Atomic Data 



1 INTRODUCTION 

Several ionization stages of carbon-like ions provide useful 
diagnostics in the mid- to far-infrared (MIR-FIR) space- 
borne observations, as well as prominently observed opti- 
cal lines from ground based instruments. Nebular plasmas 
are generally the most common sources of these lines arising 
from forbidden transitions within the levels of the ground 
configuration of C-like ions (viz. Dopita and Sutherland 
2003, Pradhan and Nahar 2011). Owing to much lower ex- 
tinction than O/UV lines, the MIR and FIR reveal spatially 
extended sources to much greater depths. Therefore, in re- 
cent years space based observatories have provided added 
impetus to the analysis of MIR and FIR lines observed by 
Spitzer, Herschel and Sofia. In addition to molecular clouds 
photoionized by hot stars into nebulae, a much wider variety 
of astronomical objects whose physical conditions may be 
explored with mid-IR [Nev] lines. These range from active 
galactic nuclei (Melendez etal. 2011, Perez-Beaupuits etal. 
2011, Stern etal. 2012, Smith etal. 2012), metal-rich Hll 
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regions (Furness etal. 2010), dusty environments of Ultra- 
Luminous-Infrared-Galaxies (ULIRGs, Houck etal. 2004, 
2005; Nagao etal. 2011), and blue compact dwarf galaxies 
(Izotov etal. 2012). 

H II regions in general, and planetary nebulae (PNe) in 
particular, are the testing ground as well as points of applica- 
tion of forbidden line ratios diagnostics with respect to tem- 
perature, density, and chemical abundances. The PNe are 
relatively bright and apparently simple objects that act as 
astrophysical laboratories to study basic atomic processes. 
However, several previous studies have discussed problems 
regarding [Nev] lines (Clegg etal. 1987, Oliva etal. 1996, 
van Hoof et al. 2000). Based on the MIR 14/24 /im line ratio, 
Rubin etal. (2001) and Rubin (2004) pointed out that for a 
number of PNe the observed values lie much lower than the- 
oretically possible emissivity ratios at standard nebular tem- 
peratures and densities. More recently, in studies of spectral 
energy distributions in active galactic nuclei (AGN) using 
high-ionization mid-IR lines, a number of observations also 
reveal anomalous 14/24 /im line ratios found to be below 
the theoretical limit (viz. Melendez etal. 2011, Tominasin 
etal. 2010, Weaver etal. 2010, Pereira-Santaella etal. 2010). 
Although the line forming regions in Seyfert galaxies may 
be dust obscured, the lower than predicted ratios are not 
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countenanced by derived mid-IR extinction (Li and Draine 
2001). It is important to address this issue in order to use 
observed [Ne v] line ratio as density diagnostics, as the emis- 
sivity ratios indicate. Moreover, a deviation from theoretical 
hmits of the otherwise clean and uncluttered spectral [Ne v] 
lines also renders the determination of mid-IR extinction 
uncertain (Melendez etal. 2011). 

The analysis of the MIR and FIR space observations 
still rests on our understanding of nebular astrophysics. For- 
bidden lines are the primary diagnostics for low-temperature 
and low-density plasmas. In principle, the atomic physics is 
relatively straightforward: electron impact excitation of low- 
lying levels followed by radiative decay via magnetic dipole 
or electric-quadrupole transitions. These are the only two 
parameters required for simple coUisional-radiative mod- 
els, usually involving no more than the few levels of the 
ground configuration. However, it is well known that the 
electron impact excitation collision strengths and the rele- 
vant transition probabilities need to be calculated to high 
accuracy in order to determine the temperature-density 
regimes in the source. For example, C-like ground config- 
uration ls^2s^2p^ gives rise to 3 LS terms and 5 fine struc- 
ture levels '^Po,i,2,^ D2,^ So- In an earlier study, Palay etal. 
(2012) showed that near-threshold resonances in the collision 
strengths of C-like 03 make significant contribution to fine 
structure transitions. Rather elaborate computations includ- 
ing relativistic and resonance effects are necessary in order 
to delineate resonances and obtain accurate rate coefficients. 
In this work, we study the transitions in the Nev ion and 
find similar eflFects. As Rubin etal. (2001) and Rubin (2004) 
have noted, different sets of collision strength data avail- 
able in literature (viz. Lennon and Burke 1991, 1994; Griffin 
and Badnell 2000) did not resolve the discrepancy in the 
14/24 ratio lying below the theoretical low density limit. 
The aim of this paper is to re-examine the Nev collision 
strengths and line ratios computed using the Breit-Pauli-R- 
Matrix (BPRM) method (Berrington etal. 1995), with par- 
ticular emphasis on the MIR lines under low-energy-density 
(LED) conditions. In order to examine the theoretical be- 
haviour, we also extend the calculations to to Tc< lOOOK, 
not considered in previous works. Based on the more accu- 
rate data computed herein, the results should enable revised 
interpretation of MIR [Ne v] lines and originating plasma en- 
vironments. 



2 THEORY AND COMPUTATIONS 

A brief theoretical description of the BPRM calculations is 
as follows (a discussion of the general methodology and ap- 
plication to atomic processes is given, for example, in Prad- 
han and Nahar 2011). 



2.1 The Breit-Pauli approximation 

The relativistic BPRM Hamiltonian is given as 

rrBP _ 

-"jv+i — 

EN+l J Y72 2Z I v-^JV+1 2 \ , rrmass , rrDar , rrso 

i=l \~^^ - —+2^j>i —f+^N+l+^N+l+-tiN+l- 

where the last three terms are relativistic corrections, re- 
spectively: 



Table 1. Levels and energies (Et) of target (core ion) Ne V 





Level 


Jt 


EtiRy) 
NIST 


EtiRy) 
SS 


1 


Is22s22p2(3p) 





0.0 


0. 


2 


Is22s22p2(3p) 


1 


0.003758 


0.0030391 


3 


Is22s22p2(3p) 


2 


0.010116 


0.011366 


4 


ls22s22p2(iD) 


2 


0.276036 


0.30391 


5 


Is22s22p2(i5) 


2 


0.582424 


0.57413 


6 


ls22s2p3(5S'°) 


2 


0.8052 


0.71604 


7 


ls22s2p3(3£)o) 


3 


1.60232 


1.62957 


8 


ls22s2p3(3Z)°) 


2 


1.60296 


1.62932 


9 


ls22s2p3(3Do) 


1 


1.60316 


1.62929 


10 


Is22s2p3(3p°) 


2 


1.89687 


1.92363 


11 


ls22s2p3(3po) 


1 


1.89687 


1.92340 


12 


ls22s2p3(3po) 





1.89719 


1.92328 


13 


ls22s2p3(lD°) 


2 


2.46556 


2.59326 


14 


ls'^2s2p^i^S°) 


1 


2.54576 


2.64956 


15 


Is22s2p3(ip°) 


1 


2.76854 


2.88988 


16 


Is22s2p4(3p) 


2 


3.76063 


3.86076 


17 


Is22s2p4(3p) 


1 


3.76778 


3.86807 


18 


Is22s2p4(3p) 





3.77085 


3.87155 


19 


lsHs2p\^D) 


2 




4.13816 


20 


ls^2s2p'^{^S) 







4.74472 



the mass correction term, iy™^''= — —^ T^^pt, 

the Darwin term, H^'"' = ^ Ei '^^(t"). (2) 

the spin — orbit interaction term, H^° — Zo? E- ^li.Si. 

Eq. (2) representes the one-body terms of the Breit in- 
teraction. In addition, another version of BPRM codes in- 
cluding the two-body terms of the Breit-interaction (Prad- 
han and Nahar 2011; W. Eissner and G. X. Chen, in prepa- 
ration) has been developed, and is employed in the present 
work. 



2.2 Ne V target representation 

The coupled channel method embodied in the BPRM ap- 
proximation is crucially dependent on the accuracy of the 
wavefunctions of the N-electron target ion included to rep- 
resent the (N+l)-electron (e -|- ion) system. We employ 
the general purpose atomic structure code SUPERSTRUC- 
TURE (Eissner et al. 1974) to compute a wavefunction ex- 
pansion including the first 20 fine structure levels. Table 1 
compares the calculated eigenenergies with experimental 
values tabulated by the U.S. National Institute of Standards 
and Technology (www.nist.gov). 

The spectroscopic levels and configurations included in 
the (e -I- Nev) wavefunction expansion are as in Table 1. 
The full configuration-interaction basis set optimized with 
SUPERSTRUCTURE is (Eissner etal. 1974, Nahar etal. 
2003): [ls2]2s^2p^ 2s2p^ 2s^2p3s, 2p'', 2s^2p3p, 2s'^2p3d, 
2s^2p4s, 2s^2p4p, 2s2p'^3s, 

2s2p'^3p, 2s2p'^3d, 2s^3s^,2s'^3p'^, 2s^i(f ,2s^As^ ,2s^Ap^ , 
2s^2,s2,p, 2s^3s4s, 2s^3p3d, 2p^3s, 2p^3p, 2p^M. Although the 
computed enegies are generally within a few percent of the 
observed values, the latter are used in the BPRM calcula- 
^ -tions to ensure precise resonance positions relative to the 
Nev target thresholds. The collision strengths were com- 
puted employing the extended BPRM codes including a full 
representation of the two-body Breit terms (Eissner and 
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Figure 1. Rydberg series of near-threshold resonances in the 
collision strength for the fine structure transition 2p^('^Pa — ^ Pi) 
corresponding to the 24 /^m FIR line. A fine mesh order 10~® 
Ryd was used to fully resolve the {^P2,^ ^2i^ So)nl resonances, 
necessary for accurate rate coefficients at low temperatures. 



Chen 2012). It is important to ensure convergence of col- 
lision strengths with respect to partial waves and energy 
resolution. Total (e 4- ion) symmetries up to (LS)j7r with J 
^ 19.5 were included in the calculations, though it was found 
that the collision strengths for all forbidden transition tran- 
sitions converged for J ^ 9.5. A very fine energy mesh of 
AE < 10~® Rydbergs was used to resolve the near-thresold 
resonances. 



2.3 Effective collision strengths 

The effective collision strengths or rate coefficients are 
computed by convolving the collision strengths over a 
Maxwellian function at a given electron temperature Teas 



r^ATe) 



f7y(e)exp(-e/fcTe)d(e/fcre), 



(3) 



where Eij is the energy difference and Q,ij is the collision 
strength for the transition i -^ j. Later we discuss LED 
plasma conditions where the effective collision strength may 
differ from the form in Eq. (3). The excitation rate coefficient 
is related to the effective collision strength as 

„ fT^ \ 8-63 X 10 ^-EjjkT^r^ ,rpx ,.-. 

qijiTe) = —17:5 — e "' Ty(re), (4) 



g^T, 



1/2 



where gi is the statistical weight of the initial level. 
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Figure 2. Collision strengths over extended energy ranges for 
fine structure transitions 2p^(^Po — ^ Pl,^ Pi —^ P2) at a) 24 /^m 
b) 14 fim respectively. 
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Figure 3. Maxwellian averaged effective collision strengths T{Te) 
for the transitions 2p2 (3 Po- 3 Pi, ^ Pi -^P2) at a) 24 ^m b) 14 fim 
(cf. Fig. 2):solid line. Dotted lines show previous results without 
relativistic effects (Lennon and Burke 1994); dash-dot and dashed 
lines show 20-level BPRM and 138-level ICFT results respectively 
(Griffin and Badnell 2000), available in the temperature range 
Te > 1000 K. The enhancements in present calculations are due 
to the fully resolved resonance structures near the threshold (cf. 
Fig. 1). 



3 RESULTS AND DISCUSSION 

Collision strengths and line ratios are presented in this sec- 
tion. A discussion of possible implications for LED plasmas 
in general, and PNe in particular, is also given. 



3.1 Ground state fine structure collision strengths 

One of the main aims of the present calculations was com- 
plete resolution of resonance structures in near-threshold 
collision strengths. Fig.[T]shows the resolved Rydberg series 
of resonances in Q{^Po —^ Pi) converging on to the higher 
-P2, D2, So levels. The calculations were carried out at an 
energy interval of 10~^ Ryd. These resonances dominate the 
low-energy-temperature behaviour of the detailed, as well as 



the averaged, collision strengths. The corresponding 24 fim 
line is the one most affected. Fig.[2]shows the results over an 
extended energy range for r2(^Po -^ Pi) and ^(^Pi -^ P2). 
Particularly noteworthy is the significant rise close to thresh- 
olds for both transitions. 

The Maxwellian averaged collision strengths are plot- 
ted in Fig. [3] (solid curve), and compared with previous 
results from Lennon and Burke (1994, dotted curve) and 
two sets of calculations by Grfiffin and Badnell (2000). The 
latter work employed two different approximations, a 20- 
level BPRM calculation (dash-dot curve) and a more ex- 
tensive intermediate-coupling fine structure transformation 
(ICFT, dashed curve) that is essentially non-relativistic (as 
in Lennon and Burke 1994) but incorporates relativistic 
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terms perturbatively via an algebraic transformation. For 
the 24 /xm transition, we find differences of up to 30% at 
the lowest temperatures with the otherwise similar previous 
20-level BPRM results, most likely attributable to improved 
resolution at very low energies in the present calculations. 
Differences no larger than 10% are found for the 14 ^m tran- 
sition between all sets of calculations. 

Present calculations are also carried out at Te< lOOOK, 
not considered in previous works. Since we have taken 
particular care to fully resolve the near-threshold collision 
strengths, the effective collision strengths should be limited 
only by the intrinsic accuracy of the BPRM calculations and 
not by numerical resolution of resonances. That is neces- 
sary in order to interpret the line ratio observations in cases 
where they do not fall in the nebular temperature range, as 
discused next. 



3.2 The 14/24 /im emissivity ratio and 
observations 

Owing to small differences in excitation energies, the fine 
structure line ratios due to transitions within the ground 
state levels ■^Po,i,2 are not highly temperature sensitive in 
the typical nebular range around 10,000K. On the other 
hand, this emssivity ratio is highly density sensitive at typ- 
ical nebular densities of 10''"® cm"''. Rubin (2004) has tab- 
ulated the 14/24 line flux ratio from a number of PNe, and 
noted that it falls below the theoretically possible value of 
approximately 1.0 in about half the observations. Fig. 0] 
shows the measured line ratios; the asterisks denote the de- 
rived densities for the PNe that fall very well on the 10,000 
K emissity ratio curve. However, the observed values from 
some of the anomalous PNe are seen as lying around the 
emssivity ratio curve at lOOOK or lower, shown as vertical 
dashed band in Fig. [J] One can readly find low temperature 
curves vs. density that could account for these ratios lying 
below the 10,000K curve. Incidentally, variations by similar 
amounts in both transitions yield very good agreement with 
ratios at Te> 10,000K with previous works (viz. Lennon and 
Burke 1994), shown in the right panel. 

Based on this work we can rule out further uncertainties 
in the collision strengths. Fig.[4]illustrates the dichotomy be- 
tween the PNe with normal density sensitivity of the 14/24 
ratio, and a fair number of those that appear to correspond 
to much lower temperatures than that which characterizes 
typical nebular photoionization equilibrium. While the equi- 
librium models may not explain these values, it is worth 
examining the behaviour of line ratios in LED plasmas in 
general. Eq. (3) implies that 



Zz77iT->oT(r) — liniE^o^iE). 



(5) 



Between 1000-10, OOOK the Mawellian factor involves ex- 
ponential factor exp(-E/kTe) ranging from exp(-158 E) to 
exp(-15.8E) in Rydberg units. Therefore T(r) depends en- 
tirely on the detailed behaviour of n{E). Furthermore, if 
we also consider that fact that low densities Ne~ 10^ cm~'^ 
may exist, a departure from the normal Maxwellian distri- 
bution under LED conditions is possible. Therefore, the rate 
coefficient would tend to follow the collision strengths and 
variations with resonance structures at very low energies 
close to exciation threshold(s). This suggestion is offered to 
re-examine non-equilibrium conditions that might prevail in 
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Figure 4. The density and temperature dependence of the 14/24 
lira line emissivity ratios: a) Solid lines represent line ratios at four 
temperatures. Asterisks denote observed line ratios from planetary 
nebulae (Rubin 2004), at Te= 10,000 K with assigned densities. 
Dotted lines represent observed line ratios that are out of range 
of typical nebular temperature-density range (Rubin 2004), ex- 
cept at much lower temperatures, b) solid curves arc present re- 
sults and the dashed lines are using the Lennon and Burke (1994) 
collision strengths. Large differences in line ratios are not found 
because of systematic differences of similar magnitude for both 
lines between previous and present rate coefficients. 



LED plasma sources, which may also not have established 
a Maxwellian distribution. The emitted line intensity would 
then depend on the actual kinetic energy and number den- 
sity of electrons at a given energy. In addition, the electron 
and ion temperatures may differ, i.e. T^/ Ti. 



3.3 [Nev] optical collision strengths and line 
ratios 

Collision strengths for forbidden optical transitions within 
ground configuration levels are shown in Fig. O "^Pi — ^ D2, 
'^P2 -^ D2,^ So -^ D2 at AA 3345.8, 3425.9, 2972.8, respec- 
tively. Resonance structures due to Nev target excitation 
thresholds given in Table 1 manifest themselves at all ener- 
gies, and need to be considered in the calculation of aver- 
aged collision strengths. Whereas the agreement with pre- 
vious works (Lennon and Burke 1994, Griffin and Badnell 
2000) is found to be satisfactory, generally 10-15% or bet- 
ter, for some of the relatively weak transitions there are 
marked differences. Fig. [S] compares the present results for 
the maxwellian averaged collision strengths for two transi- 
tions, ^Po —^ D2 and ■^Po — ^ So that do not directly cor- 
respond to observed lines but enter in the solution of rate 
equations for line emissivities. The differences are found to 
be rather large. Nevertheless, the important optical line ra- 
tios are not significantly affected, as shown in Fig. [T] The 
blended emissivity ratio (3426-1-3346) /2973 differs little from 
those derived from earlier values from Lennon and Burke 
(1994). 



3.4 Maxwrellian averaged collision strengths 

In Table 1, we present the effective collision strengths (Eq. 
3) for the 10 transitions among the ground configuration 
levels and their wavelengths. The tabulation is carried out 
at a range of temperatures typical of nebular environments. 
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Table 2. Effective Maxwellian averaged collision streiiErths 



Transition 



T(IOO) T(500) T{1000) T(5000) T(IOOOO) T(20000) T(30000) 



^Po- 


-'Pi 


14 ^m 


2.72(0) 


3.02(0) 


2.78(0) 


2.09(0) 


1.69(0) 


1.28(0) 


1.07(0) 


^P2- 


-'Po 


9 ^m 


2.70(0) 


3.00(0) 


3.16(0) 


2.56(0) 


1.93(0) 


1.34(0) 


1.04(0) 


«P2- 


-■'Pi 


14 fira 


8.48(0) 


9.18(0) 


9.55(0) 


8.02(0) 


6.26(0) 


4.49(0) 


3.60(0) 


'D2 


-^Po 


3301.3 A 


1.23(-1) 


1.20(-1) 


1.22(-1) 


1.80(-1) 


1.95(-1) 


2.09(-l) 


2.19(-1) 


'D2 


-^Pi 


3345.8 A 


3.76(-l) 


3.67(-l) 


3.69(-l) 


5.37(-l) 


5.84(-l) 


6.27(-l) 


6.58(-l) 


'D2 


-''P2 


3425.9 A 


6.61(-1) 


6.46(-l) 


6.47(-l) 


9.22(-l) 


9.97(-l) 


1.07(0) 


1.1(0) 


'So- 


-^Po 


- 


4.60(-2) 


4.82(-2) 


4.84(-2) 


3.58(-2) 


3.15(-2) 


2.98(-2) 


2.94(-2) 


'So' 


-■'Pi 


1574.8 A 


1.35(-1) 


1.42(-1) 


1.42(-1) 


1.04(-1) 


9.13(-2) 


8.73(-2) 


8.69(-2) 


'So- 


-■'P2 


1592.3 A 


2.24(-l) 


2.35(-l) 


2.35(-l) 


1.73(-1) 


1.53(-1) 


1.47(-1) 


1.45(-1) 


'So- 


-'D2 


2972.8 A 


3.51(-1) 


4.51(-1) 


4.88(-l) 


6.69(-l) 


6.26(-l) 


6.33(-l) 


6.87(-l) 
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Figure 5. Collision strengths for the forbidden optical transi- 
tions 2p2(3pi -1 D2,' P2 -^ P'2,^ So -^ D2) at 3346,3426,2973 A 
respectively 
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Figure 6. Maxwrellian averaged collision strengths T(Te) for the 
transitions 2p2(3Po -^ D2? Po -^ So): Solid Unes. Dotted lines 
show previous results without relativistic effects (Lennon and 
Burke 1994); dashed and dash-dot lines show 20-level BPRM and 
138-level ICFT results respectively (Griffin and Badnell 2000) 



including the low temperature range T ^ lOOOK not hereto- 
fore considered. 



3.5 Conclusion 

New calculations for [Nev] collision strengths are reported 
including relativistic and resonance effects at very low en- 
ergies that affect the low temperature rate coefficients not 
heretofore calculated. In particular, a precise delineation of 
near-threshold resonance structures has a significant effect 
on the mid-IR 14/24 /^m line emissitivies. However, low- 
density theoretical limits of the ratios computed at typi- 
cal nebular temperatures Te~ 10,000K are still higher than 




Figure 7. Blended [Nev] optical line ratio (3346-1- 3426)/2973 vs. 
Teat Ne= 10^ cm~^. Despite significant differences in the colli- 
sion strengths the optical emissivity ratios agree with those using 
earlier data by Lennon and Burke (1994) owing to systematic 
enhancements in the same direction. 



the anomalously low values of observations from a number 
of PNe discussed by Rubin etal. (2001) and Rubin (2004). 
Since the Maxwellian averaged collision strengths generally 
decrease with temperature, and simulate the behaviour of 
collision strengths at low energies, it is suggested that non- 
equilibrium conditions corresponding to LED plasma con- 
ditions might possibly account for the low observed 14/24 
ratios. In any case, the low temperature rate coefficients for 
Te< 10,000K are likely to be uncertain owing to resonances 
at low energies. 

The [Ne v] forbidden optical collision strengths are gen- 
erally in good agreement with previous works, and the 
prominent line ratios are not much influenced. Finally, as 
noted in the previous work on 03 lines (Palay etal. 2012), 
for higher temperatures T > 20,000K proton impact exci- 
tation of the ground state fine structure levels ^Po,i,2 needs 
to be taken into account; at lower temperatures the excita- 
tion rate coefficent due to electrons far exceeds that due to 
protons (Ryans etal. 1999). 

Collision strengths for all 190 transitions among en- 
ergy levels in Table 1 may be obtained from Ethan Palay 
(palay.5@buckeyemail.osu.edu). 
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